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Abstract 

The quantum tunneling of the magnetization vector between excited levels 
are studied theoretically in single-domain ferromagnetic nanoparticles with bi- 
axial crystal symmetry placed in an external magnetic field at an arbitrarily 
directed angle in the ZX plane. By applying the periodic instanton method 
in the spin-coherent-state path-integral representation, we calculate the tun- 
nel splittings and the tunneling rates between excited levels in the low barrier 
limit for different angle ranges of the external magnetic field (6h = tt/2, 
7t/2 <C Or <C it, and Oh = tt). The temperature dependences of the tunneling 
frequency and the decay rate are clearly shown for each case. Our results show 
that the tunnel splittings and the tunneling rates depend on the orientation 
of the external magnetic field distinctly, which provides a possible experimen- 
tal test for magnetic quantum tunneling in nanometer-scale single-domain 
ferromagnets. 
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I. INTRODUCTION 



Macroscopic quantum tunneling (MQT) and coherence (MQC) of the magnetization 
were intensively investigated both theoretically and experimentally in recent years.0 More 
recently, much attention was attracted to the spin tunneling in the single-domain ferromag- 
netic (FM) nanoparticles in the presence of an external magnetic field applied at an arbitrary 
angle. The MQT problem for FM particles with uniaxial crystal symmetry was first studied 
by Zaslavskii with the help of mapping the spin system onto a one- dimensional particle 
system.! For the same crystal symmetry, Miguel and Chudnovskyi calculated the tunneling 
rate by applying the imaginary-time path integral, and demonstrated that the angular and 
field dependences of the tunneling exponent obtained by Zaslavskii's method and by the 
path-integral method coincide precisely. Kim and Hwang performed a calculation based 
on the instanton technique for FM particles with biaxial and tetragonal crystal symmetry.0 
Kim extended the tunneling rate for biaxial crystal symmetry to a finite temperature, and 
presented the numerical results for the WKB exponent below the crossover temperature and 
their approximate formulas around the crossover temperatureJi The quantum-classical tran- 
sition of the escape rate for FM particles with uniaxial crystal symmetry in an arbitrarily 
directed field was studied by Garanin, Hidalgo and Chudnovsky with the help of mapping 
onto a particle moving in a double-well potential.! The switching field measurement was 
carried out on single-domain FM nanoparticles of Barium ferrite (BaFeCoTiO) containing 
about 10 5 — 10 6 spins.! The measured angular dependance of the crossover temperature was 
found to be in excellent agreement with the theoretical prediction,! which strongly suggests 
the MQT of magnetization in the BaFeCoTiO nanoparticles. Lii et al. studied the MQT 
and MQC of the Neel vector in single-domain antiferromagnetic (AFM) nanoparticles with 
biaxial, tetragonal, and hexagonal crystal symmetry in an arbitrarily directed field.! 

It is noted that the previous results of MQT of the magnetization vector at excited 
levels in an arbitrarily directed field were obtained by numerically solving the equation of 
motion satisfied by the least trajectory.! The purpose of this paper is to present an analytical 
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investigation of the quantum tunneling at excited levels in the biaxial FM particles in an 
arbitrarily directed field, based on the periodic instanton method.1'0 Both the nonvacuum 
(or thermal) instanton and bounce solution, the WKB exponents and the preexponential 
factors are evaluated exactly for different angle ranges of the magnetic field (Oh = tt/2, 
7r/2 + 0(e 3 / 2 ) < H < 7T — 0(e 3 / 2 ), and Oh = tt). Our results show that the distinct angular 
dependence, together with the dependence of the WKB tunneling rate on the strength of 
the external magnetic field, may provide an independent experimental test for the magnetic 
tunneling at excited levels in single-domain FM nanoparticles. 

This paper is structured in the following way. In Sec. II, we review briefly some basic 
ideas of MQT and MQC in FM particles. And we discuss the fundamentals concerning the 
computation of level splittings and tunneling rates of excited states in the double-well-like 
potential. In Sees. Ill, we study the spin tunneling at excited levels for FM particles with 
biaxial crystal symmetry in the presence of an external magnetic field applied in the ZX 
plane with a range of angles n/2 < H < vr. The conclusions are presented in Sec. V. 

II. MQT AND MQC OF THE MAGNETIZATION VECTOR IN FM PARTICLES 

The system of interest is a nanometer-scale single-domain ferromagnet at a temperature 
well below its anisotropy gap. For such a FM particle, the tunnel splitting for MQC or the 
tunneling rate for MQT is determined by the imaginary-time transition amplitude from an 
initial state \i) to a final state |/) as 

Ufi = (/I e~ HT \i) = J DQexp (S E ) , (1) 

where Se is the Euclidean action and DQ is the measurement of the path integral. In the 
spin-coherent-state representation, the Euclidean action can be expressed as 



M d<p\ .Mo (d<t>\ „ w 
i — — — % — — cos 9 + E 
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(2) 



where V is the volume of the FM particle and 7 is the gyromagnetic ratio. Mq 



hjS/V, where S is the total spin of FM particles. It is noted that the first two terms 



in Eq. (2) define the topological Berry or Wess-Zumino, Chern-Simons term which arises 
from the nonorthogonality of spin coherent states. The Wess-Zumino term has a simple 
topological interpretation. For a closed path, this term equals —iS times the area swept 
out on the unit sphere between the path and the north pole. The first term in Eq. (2) is 
a total imaginary-time derivative, which has no effect on the classical equations of motion, 
but it is crucial for the spin-parity effect s.BEiHll However, for the closed instanton or bounce 
trajectory described in this paper (as shown in the following), this time derivative gives a 
zero contribution to the path integral, and therefore can be omitted. 

In the semiclassical limit, the dominant contribution to the transition amplitude comes 
from finite action solution (instanton or bounce) of the classical equation of motion. The 
instanton's contribution to the tunneling rate T or the tunnel splitting A (not including the 
topological Wess-Zumino phase) is given by0 



where u> p is the oscillation frequency in the well, S c i is the classical action, and the prefactor 
A originates from the quantum fluctuations about the classical path. It is noted that Eq. (7) 
is based on quantum tunneling at the level of ground state, and the temperature dependence 
of the tunneling frequency is not taken into account. However, the instanton technique is 
suitable only for the evaluation of the tunneling rate or the tunnel splitting at the vacuum 
level, since the usual (vacuum) instantons satisfy the vacuum boundary conditions. Recently, 
Liang et a/.&0 developed new types of pseudoparticle configurations which satisfy periodic 
boundary condition (i.e., periodic instantons or nonvacuum instantons). They found that 
the tunneling effect indeed increases exponentially with energy in the low-energy region. 

For a particle moving in a double- well-like potential U (x), the level splittings of degen- 
erate excited levels or the imaginary parts of the metastable levels at an energy E > are 
given by the following formula in the WKB approximation,!!^!^ 




(7) 




(8) 
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and the imaginary-time action is 



/ 

Jx 



■x 2 (E) 



dx^/U (x) - E, 



(9) 



S (E) = 2V2m 



xi(E) 



where X\^ (E) are the turning points for the particle oscillating inside the inverted potential 
—U (x). uj (E) = 2n/t (E) is the energy- dependent frequency, and t (E) is the period of the 
real-time oscillation in the potential well, 



where 23,4 (E) are the turning points for the particle oscillating inside the potential U (x) . 
Recently, the crossover from quantum to classical behavior and the associated phase transi- 



In this section, we study the tunneling behaviors of the magnetization vector in single- 
domain FM nanoparticle which has the biaxial crystal symmetry, with ±z being the easy 
axes in the absence of an external magnetic field. The magnetic field is applied in the ZX 
plane, at an angle in the range of n/2 < 9h < tt. Then the magnetocrystalline anisotropy 
energy E (9, 0) can be written as 

E (9, 0) = K x sin 2 6 + K 2 sin 2 9 sin 2 - M H X sin 6 cos <p - M H Z cos 9 + E , (11) 

where K\ and K 2 are the longitudinal and the transverse anisotropy coefficients, respectively, 
and Eq is a constant which makes E (9, 0) zero at the initial orientation. As the external 
magnetic field is applied in the ZX plane, H x = H sin 9h and H z = H cos 9h, where H is 
the magnitude of the field and 9h is the angle between the magnetic field and the z axis. 
By introducing the dimensionless parameters as 




(10) 



tion were studied extensively in nanospin system: 



JlJiZHi] anc [ other systems.ill 



III. MQC AND MQT FOR BIAXIAL CRYSTAL SYMMETRY 



K 2 = K 2 /2K lt H x = H x /H , H z = H z /H ( 



(12) 



the E (9, 0) term of Eq. (11) can be rewritten as 
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E (0, 0) = - sin 2 + K 2 sin 2 sin 2 - H x sin cos - if ^ cos + £ , 



(13) 



where E (0, 0) = 2K\E (0, 0), and if = 2Ki/M . At finite magnetic field, the plane given 
by = is the easy plane, on which E (9, 0) reduces to 

1 



£(0,0 = 0) = -sin 2 0-# cos (9 - 9 H ) + E . 
2 



(14) 



We denote 0o to be the initial angle and 9 C the critical angle at which the energy barrier 
vanishes when the external magnetic field is close to the critical value H c {Oh) (to be calcu- 
lated in the following). Then, O satisfies dE {0, = 0) jdO _ = 0, C and H c satisfy both 



dE {0, = 0) jdO 



=0 and 

,H=H C 



d 2 E {0, = 0) jd0 r < 



=0, which leads to 



,H=H 



-sm(20 o ) + Hsm(0 o -0 H ) = 0, 
^sm(20 c ) + H c sm(0 c -0 H ) = 0, 
cos {20 c ) + H c cos {0 C - H ) = 0. 



(15a) 
(15b) 
(15c) 



After some algebra, the dimensionless critical field H c {Oh) and the critical angle C are found 
to be 



sin (20 c ) 



(sin0 H ) 2/3 + |cos0^| 2/3 
2|cot H \ 1/3 



-3/2 



(16a) 
(16b) 



1 + |cot^| 2/3 ' 

Now we consider the limiting case that the external magnetic field is slightly lower than 
the critical field, i.e., e = 1 — H/H C <C 1. At this practically interesting situation, the barrier 
height is low and the width is narrow, and therefore the tunneling rate in MQT or the tunnel 
splitting in MQC is large. Introducing rj = C — {\rj\ <C 1 in the limit of e 1), expanding 



dE {9 
and 



0)/d0 



about 6 a and using the relations dE (0, = 0) jdO 



M=h c 



d 2 E (0, = 0) /d0 2 ___ = 0, Eq. (15a) becomes 



,H=H 
3 



sin (20 c ) (e - -rf^j - rj cos (20 c ) (2e - rf) = 0. 



(17) 



Then the potential energy E (0, 0) reduces to the following equation in the limit of small e, 



E (5, 0) = K 2 sin 2 sin 2 (0 O + 6) + H x sin (0 O + S) (1 - cos 0) + ^ (5) , 



(18) 



where 5 = 9 — 9 (\S\ <C 1 in the limit of e 1), and E t (5) is a function of only 5 given by 

(19) 



E 1 (5) = j sin (29 c ) (35 2 r] - 5 3 ) + ]- cos (20 c ) 



^ ( < - ^ 2 ) + <^ - ^ 



In the following, we will investigate the tunneling behaviors of the magnetization vector 
at excited levels in single-domain FM nanoparticles with biaxial crystal symmetry at three 
different angle ranges of the external magnetic field as 9h = tt/2, tt/2 + O (e 3 ^ 2 ) < &h < 
7r — O ( e3 ^ 2 )> and ®h = 7r, respectively 



A. 6 H = tt/2 

For 9h = 7r/2, we have C = 7r/2 from Eq. (16b) and rj = \f2e from Eq. (17). Eqs (18) 
and (19) show that is very small for the full range of angles n/2 < 9 H < n for FM particles 
with biaxial crystal symmetry. Performing the Gaussian integration over 0, we can map the 
spin system onto a particle moving problem in the one-dimensional potential well. Now the 
imaginary-time transition amplitude Eqs. (1) and (2) becomes 



U/i = J d5 exp (S E [5]) , 
= J d5 exp \—jdr 



with 



1 fdSY TT . £ . 
2™ U +U{6) 



hs 2 



(20) 



m 



2V[K 2 + K 1 (l-e)]' 



and 



»(fi = ff(i-^'. 



(21) 



The problem is one of MQC, where the magnetization vector resonates coherently between 
the energetically degenerate easy directions at 5 = and 5 = 2y / 2e separated by a classically 
impenetrable barrier at 5 — y^2e. 
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Now we apply the periodic instanton methodEHlj to evaluate the level splittings of excited 
states. The periodic (or thermal) instanton configuration 5 P which minimizes the Euclidean 
action in Eq. (20) satisfies the equation of motion 

2 



Hf) - u{6 > ] - e ' 



(22) 



where E > is a constant of integration, which can be viewed as the classical energy of the 
pseudoparticle configuration. Then the kink-solution is0 



5 P = v2e + V2e — asn (&it, k) 



(23) 



where a = 2JhE/KiV, and uui = J KiV/2fmi\ / 2e + a. sn(u)iT,k) is the Jacobian elliptic 



sine function of modulus k = J (2e — a) / (2e + a). In the low energy limit, i.e., E — > 0, 



k — > 1, sn(w, 1) — > tanhw, we have 



5 p = v 2e [1 + tanh (^it)] 



(24) 



which is exactly the vacuum instanton solution derived in Ref. 4, where uj\ 



V2 {K x V/hS) e 1 / 2 ^ - e + K 2 /K x . 

The Euclidean action of the periodic instanton configuration Eq. (23) over the domain 
(—(3, (3) is found to be 

y 

2"" V rfr 



/3 



im[%] + 17 (dip) 



W + 



(25) 



with 



W 



] K x Vm 



(26 



.3/2 



,/2 



2 - k' 2 



K{k) 



(26) 



3\/2V n jl-k»/2l 

where fc /2 = 1 — k 2 . K (k) and .E (k) are the complete elliptic integral of the first and second 
kind, respectively. The general formula Eq. (8) gives the tunnel splittings of excited levels 

as 



AE=^ exp (-W) , 



7T 



(27) 



where W is shown in Eq. (26), and 00(E) = 2n/t(E) is the energy-dependent frequency. 
For this case, the period t (E) is found to be 

rS 2 



t(E) 



'2m 



i 



d5 



2hm 



1 



a 



-K (k') , 



(28) 



'Si ^E-U (5) V^i^V^T 

where 5\ = V2e + \/2t — a, and 82 = \/2e + \/2e + a. Now we discuss the low energy limit 
where E is much less than the barrier height. In this case, k' A = 4fiE / K x Ve 2 <C 1, so we 
can perform the expansions of K (k) and E (k) in Eq. (26) to include terms like k' 4 and 
k' 4 In {4/k'), 



E(k) =1 + - 



h 4 



k u + — 
16 



ff(*)=h.(±) + i[l»(± 



k'J 12 
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6 



k ■■■ . 



(29) 



With the help of small oscillator approximation for energy near the bottom of the potential 



well, E n = (n + 1/2) Q u Q ± = ^U" (5 = V^) /m = 2^K l Ve/hm, Eq. (26) is expanded as 



1 / h f 1 

64e 3 / 2 V K x Vm V + 2 



(30) 



Then the general formula Eq. (8) gives the low-lying energy shift of n-th excited states for 
FM particles with biaxial crystal symmetry in the presence of an external magnetic field 
applied perpendicular to the anisotropy axis {6 H = rr/2) as 

2ll/2 e 3/2^ \ n+1/2 



O / OH/2,3/2 q \ • 



x cxp 



2 5/2 



c 3/2 



3 VI -e + A J ' 
where A = K 2 / K 1 . 

When n = 0, the energy shift of the ground state is 

915/4 / 95/2 

hAE = l — (K.V) (l-e + A) 1/4 e 5 /4 >s -i/2 exp - 

V n \ 

Then Eq. (31) can be written as 

hAE n = ^ (hAE ) , 



(31) 



,3/2 



3 VI - e + A 



(32) 



(33) 



where 

2 11 /2 e 3/2 iS 

* = VT^7T=x- (34) 

To see the temperature dependence we take the Boltzmann average of the tunneling fre- 
quency / = AAE at temperature T, 

/ ( T ) = V- E 4A ^ ex P (-*>E n P) , (35) 

where Z Q = I] n exp (—hE n {3) is the partition function with the harmonic oscillator approxi- 
mated eigenvalues E n — (n + 1/2) f^. The final result of the tunneling frequency at a finite 
temperature T is found to be 

/ (T) = 4AE (l - e-^ l/3 ) exp ( qi e~ m ^) , (36) 

where AE and qi are shown in Eqs. (32) abd (34). 

B. vr/2 + O (e 3 / 2 ) < H < tt - O (e 3 / 2 ) 

For 7r/2 + O (e 3/2 ) < 6 H < n -O (e 3/2 ) , the critical angle 9 C is in the range of O (e 3/2 ) < 
# c < 7r/2 — (e 3 ^ 2 )) and ?7 ~ \j2e/?>. Now the problem can be mapped onto a problem of 
one-dimensional motion by integrating out 0, and for this case the effective mass m and the 
potential U (5) in Eq. (20) are found to be 

hs 2 

m = 

2K{V 

and 



l+|cot4| 2/3 + A 



U{5) =^fsm26 c (V6~e5 2 - 5 3 ) 



2h 

3U q 2 (l-^q), (37) 



where q = 35/2VQe, and U = (2 5/2 /3 3/2 ) (K-^Ve^ 2 /?^ sin2# c . The problem becomes one 
of MQT, where the magnetization vector escapes from the metastable state at 5 = 0, = 
through the barrier by quantum tunneling. 



Now the periodic bounce configuration with an energy E > is found to be 



5 P = -VQe a — (a — b) sn 2 (u; 2 r, k) 



(38) 



where 



" 2 = 2V4x3V4 V 



(39) 



a (E) > b (E) > c (E) denote three roots of the cubic equation 



<7 3 - -q 2 + — 



= 0. 



(40) 



sn(cj 2 T, k) is the Jacobian elliptic sine function of modulus k = ^J(a — b) / (a — c). In the 
t, i.e., .E — > 0, fc — > 1, sn(w, 1) — > tanhw, a — > 3/2, 6 — * 0, we have 



low energy limit, 



COSH (U^Tj 



(41) 



where 



" 2 " 2VI 



| cot 
1 + I cot H | 

bounce solution obtained in Ref. 4. 



Eq. (41) agrees well with the vacuum 

The classical action of the periodic bounce configuration Eq. (38) is 



S p — 



-13 



dr 



W + 2E/3, 



(42) 



with 



o9/2 

W = ^^^ndJo(a - cf 2 [2 (V - k 2 + l) E (k) - (l - k 2 ) (2 - k 2 ) K (k)] . (43) 



The period t (E) of this case is found to be 

d5 



t(E) 



'2m 



= 4, 



2 em 



-K(k'), 



(44) 



yjE-U (S) V 3C/ o(a-c) 

where k' 2 = 1 — A; 2 . Then the general formula Eq. (8) gives the imaginary parts of the 
metastable energy levels as 
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lmE = ^— Uxp(-WO, 



7T 



(45) 



where to (E) = 2n/t (E), and W is shown in Eq. (43). 

Here we discuss the low energy limit of the imaginary part of the metastable en- 
ergy levels. For this case, E n = (n + l/2)fi 2 , ^2 = ^U" (5 = 0) ~]m = (3/2) ^U /me, 
a « (3/2) (1 - k' 2 /A), b w (3A; ,2 /4)(l + 3A; /2 /4), c « - (3A; ,2 /4) (1 + k> 2 /4), and k' 4 = 
16E/27U < 1. Therefore, Eqs (43) and (45) reduce to 



w = 217/4 _ x3 ' /4 S ^ 



+ ( n + - ) In 



I cot 0# 



1/6 



^l-e + A(l + |cot^| 2/3 ) 



225/4 


x ?, ll l 4 Se>l 4 \cot6 H 


1/6 




/l - e + A (l + |cot0 H 


2/3^ 



(46) 



and 



,1/6 



MmE n = 3V ', X | 3/4 e^S- 1 (K.V) - ^^L_ Jl - 6 + A (l + |cot ^| 2/3 ) 



/ 



x 



2 25/4 x 3 ll/4g e 5/4| cot g H |l/6 
- 6 + A (l + I COt 0ffj 



1 + \cot 9 H \ 
\ n+1/2 



2/3 



x exp 



)17/4 



31/4 



3Ot0 



1 1/6 



i* I 



V 



1 - e + A (l + \cotO H \ 



2/3 



(47) 



For vacuum bounce case n = 0, we have 
313/9 x 2 3i/s 



ft Im £ = 



x exp 



{K X V)JI*S-W |cot ^ |1/6 



2 17/4 x 3I/4 



5e 5/4 



1 + |cot0 H | 2/3 

|cot0 H | 1/6 



l-e + A(l + |cot^| 2/3 ) 
\ 



1/4 



l-e + A(l + |cot^| 2/3 ) J 



(48) 



Then Eq. (47) can be written as 

hlmE n 

where 



n! 



(Mm£ ), 



(49) 
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2 25/4 x3 U/4S e 5/4| cot ^|l/6 

Q2 = I - (50) 

^/l-e + A(l + |cot^| 2/3 ) 

At finite temperature T the decay rate T = 2 Im E can be easily found by averaging over 
the Boltzmann distribution 

F CO = J- E Im^exp (-£ n /3) , (51) 

where Zo = J2 n ex P ( — hE n /3) is the partition function with the harmonic oscillator ap- 
proximated eigenvalues E n = (n + 1/2) f2 2 - The final result of the decay rate at a finite 
temperature T is found to be 

r (T) = 2 Im E (l - e - m ^) exp (g 2 e-^ 2/3 ) , (52) 

where lmE Q and g 2 are shown in Eqs. (48) abd (50). 

In Fig. 1 we plot the temperature dependence of the tunneling rate for the typical values 
of parameters for nanometer-scale single-domain ferromagnets: S = 5000, e = 1 — H/H c = 
2 x 10~ 3 , A = K 2 /Ki = 10, and 9h = 37r/4. From Fig. 1 we easily see the crossover 
from purely quantum tunneling to thermally assisted quantum tunneling. The temperature 
Tq '* characterizing the crossover from quantum to thermal regimes can be estimated as 
ksT^ = AU/Sq, where AU is the barrier height, and So is the WKB exponent of the 
ground-state tunneling. It is shown that in the cubic potential (q 2 — q 3 ), the usual second- 
order phase transition from the thermal to the quantum regimes occurs as the temperature 
is lowered.0 The second-order phase transition temperature is given by ksT^ = htUb/271, 



where U}, = y\U" /m is the frequency of small oscillations near the bottom of the 
inverted potential —U (x), and Xf, corresponds to the bottom of the inverted potential. For 
the present case, 5b = 2v^6e/3, 

I cot 6} 



hu b = 2 5 / 4 x 3 1 / 4 {K X V) g-i e i/4 |CUL,7h|1/ ' /l - e + A (l + |cot ^| 2/3 

1 + I cot 8h 1 ' " " 



Sq _ 2 17 / 4 x 3^ c5/4 |cot^| 1/6 



1 - e + A [1 + \cot9 H 
12 



1 2/3 



and 



3 3 ' 2 ' ' i + |cote„| 2/3 ' 

Then it is easy to obtain that 

k B T^ = X (K lV ) S-h^ |COt ^' 1/6 2/3 A /l - e + A (l + |cot ^| 2 / 3 ) , 

i + |cot^r /3 v V ' 

and k B T^ ] = (5tt/18) k B T^ 2) « 0.87A; B T (2) . 

C. Off = 7T 

In case of = 7r, we have # c = and r\ = 0. Working out the integration over 0, 
the spin tunneling problem is mapped onto the problem of a particle with effective mass 
m = hS 2 /2K 2 V moving in the one dimensional potential well U (5) = (KiV/h) (e5 2 — <5 4 /4). 
Now the problem is one of MQT, and the nonvacuum bounce at a given energy E > is 
found to be 

5 p = VTe(l + Jl--^-\ dn (u 3 r,h), (53) 



KtVe 2 



where 



1/2 



hm \ V KiVe 2 



2 

I 



1 + 

In the loe energy limit, i.e., E — > 0, k — > 1, dn(w, 1) — > 1/ cosh-u, we have 

p cosh (c<J3r) ' 



(54) 



where 0)3 = 2\[K^K^eV jTiS . Eq. (54) is in good agreement with the vacuum bounce solution 
derived in Ref. 4. 

The classical action of the nonvacuum bounce Eq. (53) is 
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rP 

S„ = dr 
J-p 



1 fdS p 
-m — — 

2 I dr 



W + 2Ep, 



(55) 



with 



A HE 



UJ 3 



(2 - k 2 ) E (k) - 2k' 2 K (k) 



(56) 



where k — 1 — k . Then the imaginary parts of the metastable energy levels are 

lmE= ^— Uxp(-W0, 

7T 

where cu (E) = 2n/t (E), and the period t (E) for this case is found to be 

1 



(57) 



t(E)=A ] 



hm 



HE 



K(k') 



(58) 



K ± Ve i + ^ 

Now we consider the low energy limit of the imaginary part of the metastable energy 



level. For this case, E n = (n + l/2)ft 3 , fi 3 = \jU" (5 = 0) ~]m = ^K^e/tim, k 
(l/2 3 / 2 e 3 / 2 ) y/h/KiVm (n + 1/2) < 1, then 

'32Se 3 / 2 ^K 1 /K. 



n+1/2 



(59) 



and 



,/9 / ^3/2 o\ n+1 /2 



exp — - 



In the case of n — 0, the imaginary part of the metastable ground state reduces to 
Then Eq. (60) can be written as 



(60) 



(61) 



n! 



(62) 



where 



?3 



32e 3 / 2 ff 
14 



(63) 



And the final result of the decay rate at finite temperature T is found to be 

r (T) = 2 Im E (l - e - nnal3 ) exp (q^- mzP ) . (64) 

The temperature dependence of the decay rate is shown in Fig. 2. It can be shown 
that the double-well potential (q 2 — q 4 ) yields the second-order phase transition from the 
thermal to the quantum regimes as the temperature is lowered. For this case, the position 
of the energy barrier is 5b = V2e, the frequency of small oscillations near the bottom of the 
inverted potential is houj, = 2 3 / 2 (y 'KiK 2 V^j e 1 l 2 S~ 1 , the WKB exponent of the ground-state 
tunneling is So = (8/3\/Xj Se 3 / 2 , and the height of barrier is hAU = (K{V) e 2 . Therefore, 
k B T^ ] = (y/2/ir) (y/KJCJf) e l ' 2 S~\ and k B T^ ] = (3tt/8v^) k B T^ 2) « 0.83£;bT (2) . 

IV. CONCLUSIONS 

In summary we have investigated the quantum tunneling of the magnetization vector 
between excited levels in single-domain FM nanoparticles with biaxial crystal symmetry 
in the presence of an external magnetic field at arbitrarily directed angle. By applying 
the periodic instanton method in the spin-coherent-state path-integral representation, we 
obtain the analytic formulas for the tunnel splitting between degenerate excited levels in 
MQC and the imaginary parts of the metastable energy levels in MQT of the magnetization 
vector in the low barrier limit for the external magnetic field perpendicular to the easy 
axis (9h =7r/2), for the field antiparallel to the initial easy axis (9h = 7r), and for the 
field at an angle between these two orientations (tt/2 + O (e 3 ^ 2 ) < H < it — O (e 3 / 2 )). The 
temperature dependences of the tunneling frequency and the decay rate are clearly shown 
for each case. One important conclusion is that the tunneling rate and the tunnel splitting 
at excited levels depend on the orientation of the external magnetic field distinctly. Even 
a small misalignment of the field with Oh = vr/2 and n orientations can completely change 
the results of the tunneling rates. Another interesting conclusion concerns the field strength 
dependence of the WKB exponent in the tunnel splitting or the tunneling rate. It is found 
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that in a wide range of angles, the e {= 1 — H/H c j dependence of the WKB exponent is 
given by e 5 / 4 , not e 3//2 for 9 H = n/2, and 9 H = it. As a result, we conclude that both the 
orientation and the strength of the external magnetic field are the controllable parameters for 
the experimental test of the phenomena of macroscopic quantum tunneling and coherence 
of the magnetization vector between excited levels in single-domain FM nanoparticles at 
sufficiently low temperatures. The theoretical calculations performed in this paper can be 
extended to the FM particles with a much more complex structure of magnetocrystalline 
anisotropy energy, such as trigonal, tetragonal, and hexagonal crystal symmetries. Work 
along this line is still in progress. We hope that the theoretical results presented in this paper 
may stimulate more experiments whose aim is observing macroscopic quantum tunneling and 
coherence in nanometer-scale single-domain ferromagnets. 
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Figure Captions: 

Fig. 1 The temperature dependence of the relative decay rate T (T) /Y (T = OK) for FM 
particles in a magnetic field with a range of angles ir/2 + O ( e3 ^ 2 ) < 6h < K — O ( e3//2 )- 
Here, S = 5000, e = 1 - H/H c = 2 x 10" 3 , A = K 2 /K x = 10, and 9 H = 3tt/4. 

Fig. 2 The temperature dependence of the relative decay rate Y (T) /Y (T = OK) for FM 
particles in a magnetic field along ir. Here, S = 5000, e = 1 — H/H c = 0.03, and A = 
K 2 jK x = 10. 
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